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O-i! Abstract 
D , 

' We discuss the spectra of multi- flavor massless QCD2. An approx- 

imation in which the Hilbert space is truncated to two currents states 
^ , is used. We write down a 't Hooft like equation for the wave func- 

^ I tion of the two currents states. We solve this equation for the lowest 

massive state and find an excellent agreement with the DLCQ results. 
In addition, the 't Hooft model and the large Nf limit spectra are 
re-derived by using a description in terms of currents. 



1 Introduction 



Two dimensional quantum chromodynamics (QCDg) is a useful toy model 
for the real world QCD. The large Nc model was solved by 't Hooft[|] and 
shows confinement of quarks with an approximately linear Regge trajectory 
of states. Other issues, such as the baryonic spectrum at strong coupling 
and questions of screening versus confinement]^, || can also be addressed in 
this framework. 

The two-dimensional model with fermions in the adjoint representation 
is also interesting and attracted a lot of attention in recent years ||, ||, 0, ||, 
^, [TD|, |TT], |r2|, In particular it was shown in that the adjoint fermions 
model is equivalent to QCD2 model with level Nf = for the massive 
part of the spectrum, in the case of massless fermions. Another attempt to 
address the adjoint fermions model, was by using the currents as building 
blocks of the spectrum [§. 

The idea of the present work is to study the spectrum of QCD2 at arbi- 
trary level Nf using states built from two currents, for the case of massless 
fermions. 

Our goal in this work is two folded: (i) To derive a 't Hooft like equation 
for the wave function of the "currentball" states for QCDg models at arbitrary 
level Nf. The equation should interpolate between the description of a single 
flavor model ('t Hooft model), the adjoint fermions model Nf = Nc and the 
large Nf model, (ii) To solve the equation for the lowest massive state. 

Whereas the 't Hooft model Nf = 1 is exactly solvable, the multi-flavor 
case Nf > 1 is not solvable model even in the Veneziano limit when both N^ 
and Nf are taken to infinity (with a fixed ratio), since pair creation and anni- 
hilation is not suppressed. In the present work , we use an approximation in 
which we restrict ourselves to two currents states. We cannot justify a-priori 
such an approximation for arbitrary level. However the numerical solutions 
for the lowest massive state admit a very close resemblance to the DLCQ 
results where such a truncation was not used |]TT], |T3|. A justification can 
be given for ^ very small or very large. 

The obtained equation suggests the following picture: the underlying 
degrees of freedom in the problem are interacting "gluons" with mass 
Actually, these are really quanta of the the color currents. As it is well known, 
there are no independent gluon degrees of freedom in two dimensions. 
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The organization of the manuscript is as follows. In section 2 the bosoniza- 
tion of multi-flavor QCD is reviewed. In section 3 we state the problem of 
diagonalizing the mass operator in terms of currents and arrive at a 't Hooft 
like equation for the currents wavefunction. Section 4 is devoted to a solution 
of the differential equation and for a discussion in some specific cases, such 
as the 't Hooft limit, the Nf S> case and the adjoint fermions [Nf = N^) 
case. Section 4 is a summary and a discussion. 



2 Massless QCD2 and Bosonization 

Massless multi-flavor QCDg with fermions in the fundamental representation 
of SU{Nc) is described by the following action 

S = Jd'x tr (-^i^. + m (1) 

where = "^i^, i = 1 . . . N^, a = 1 . . . Nf. 

It is natural to bosonize this system, since bosonization in the SU{Nc) x 
SU{Nf) X Ub{^) scheme decouples color and flavor degrees of freedom (in 
the massless case). The bosonized form of the action of this theory is given 

byl 



Sbosonized (2) 



NfSwzwih) + NcSwzwig) + J ]^d, 

J d^x tr [ih^d+hA^ + ihd^h^A+ + A+hA^h^ - A+A^ 



d^x tr -^F.^F^"' 

2g2 



Nf_ 
2tt 



where h G SU{Nc), g G SU{Nf), (p is the baryon number and Swzw stands 
for the Wess-Zumino-Witten action, which for complex fermions reads 

Swzwig) = ^ I tr {d^,gd^'g'^) + 

-I- [ d'ye^^' tr {g-%g){g-'d,g){g~'dkg), 
\2tx jb 

Since we are interested in the massive spectrum of the theory and the fla- 
vor degrees of freedom are entirely decoupled from the system and they are 
massless, we can put aside the g and cf) fields (There is a residual interaction 
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of the zero modes of the g, h and fields, but it is not important to our 
discussion^). 

Upon choosing the hght cone gauge A_ = and integrating we arrive 
to the following action 

S = NfSwzwih) - \^ J tr (3) 

where J+ = ^-^hd^h\ In terms of J = ^ the light-cone momentum 

operators take the following simple form 

= nXnj I ' -^"^^"'"^^ = 0)r(x",x+ = 0) :, (4) 
namely, the Sugawara form, and 

r 1 
p- = -_ J dx- : J"(x-, x+ = 0)^ J"(a;-, x+ = 0) : . (5) 

Our task will be to solve the eigenvalue equation 

2P+P- 1^) = M'^ (6) 

We write and P~ in terms of the Fourier transform of J{x~) defined by 
J{p^) = / '^^^ J(x~,x^ = 0). Normal ordering in the expression of 
P"*" and P^ are naturally with respect to p, where p < denotes a creation 
operator. To simplify the notation we write from now on p instead of In 
terms of these variables the momenta generators are 



= w-Tw / dpr{~p)r{p) (7) 



o2 

p^ 

Recall that the light-cone currents J'^(p) obey a level iVy, SU{Nc) affine Lie 



e 1 

P- = - / dp-J^(-p)J'^(p) (8) 

hat the light-c 

algebra 



[rip), j\p')] = p 5'^'5{p + p') + ir'^r{p + p') (9) 

We can now construct the Hilbert space. The vacuum |0,P) is defined by 
the annihilation property: 

Vp > 0, J{p) |0, P) = (10) 
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Where R is an "allowed" representations depending on the level. Therefore, 
a physical state in Hilbert space is = tr J{—pi) . . . J{—pn) |0, R). Note 
that this basis is not orthogonal. 

3 't Hooft like equation for the two currents 
wave- function 

Let us restrict ourselves to the 2-currents sector of the Hilbert space 

1^) = ^ /' dk mJ\-k)r{k - 1) |0), (11) 

namely to states which are color singlets of two currents with total = 1 
momentum and a distribution of P~ momentum Note that $ is a 

symmetric function 

= $(1 - fc). (12) 

Our task now is to find the eigenvalue (Schrodinger) equation for the 
wavefunction $(/c). Let us act with the "Hamiltonian" P~ on the state |$). 
The commutator of P~ with a current J^{—k) yields the following result 

[r^r{-p)r{p),J''{-k)] = 

Jo p^ 

(13) 

JO 

The above expression (0) contains 3 terms on the R.H.S. The first term 
contains a single creation operator. The second term contains an annihilation 
current and therefore should be again commuted with J^{k — 1). The third 
term contains two creation currents and it would lead to a 3-currents state. 
Namely, the affine Lie algebra created a higher state. This is a manifestation 
of the fact that pair creation is, generically, not suppressed in multi-flavor 
QCD2, as expected in general in QFT. 
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Note that while deriving eq.(]l3|) we get an "infinite" contribution Nc\J^{—k). 
This contribution will be canceled by a counter contribution which comes 
from the regime p ~ in the first integral on the R.H.S. of ([TB|), as below. 

The commutator of the second term in the R.H.S of (|13D with J^{k — 1) 
yields 

[ /f dp (^1 - .rV'^(-p) J^(p - kl J\k - 1)] = (14) 

r '^^ " (^^) ^"^"^-^^^"^^ -i)-r{v~ k)r{k - p - 1)). 

Our results can be summarized by the following set of equations 
1$) = 

-i- /' dk mJ''{~k).r{k - 1) |0) + (15) 

— dkdpdi 5{k + p + i-i)^{k,p,i)ir^'r{-k)j\-p)r{-i) |o) 

{NcNf)2 Jo 
with 

t: \ J 

and 

/ /I 1 \ 2N, 



-^^ r ''ii^ - L ''w^ + - (1^) lo ^^^i 

Ignoring the 3-currents term (see below), we get that $(A;) obeys the following 
eigenvalue equation 

i^l»WHiv,-iv.)(i + ^)*W 
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For general Nc and Nf discarding the 3-currents term is unjustified. However, 
since the length of is |\E'(/c,p, /)|~ e'^{NcNf)^ , in the limit of large Nc with 
fixed e^Nc and fixed Nf, or large Nj with fixed e'^Nf and fixed Nc, the 3- 
currents contribution is indeed negligible, as compared with the 2-currents 
term, the later being of order 1. Note also that while deriving eq.(|T7|) we 
assumed that Jq dp = 0. We will justify this assumption in the following. 
Another remark is that the first integral in (0) should be calculated as a 
principal value integral. The divergent part of this integral (arising from the 
regime p ^ k) cancels the previously mentioned infinity. 

In order to make contact with the ordinary 't Hooft equation, it is useful 
to integrate equation ( [T7| ) with respect to k and rewrite the equation in terms 

The derivation goes as follows. First, integrating eq.([T^) we get (p{k) = 
— k) + const. Taking ip{l) = we get 

^{k) = -^{l-k). (19) 

Now </?(!) = implies Jq dk^{k) = 0, which was our assumption above. 
Then, differentiating (|TB|) we do get (P7|), and by (|T^ we also get that there 
is no extra integration constant. 

We would like to comment on the issue of the Hermiticity of the 'Hamil- 
tonian' M^. Naively, it seems that is not Hermitian with respect to the 
scalar product < ijj\(p >= dkijj*{k)(p{k), since the Hermitian conjugate of 

However, as we shall see in the next section, the numerical solution yields real 
eigenvalues and eigenfunctions. Therefore, at least on the subspace which 
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is spanned by the eigenf unctions, namely real functions which are zero at 
A; = 0, 1 and anti-symmetric with respect to k = ^, the operator is 
Hermitian. Note that is "more regular" than (|TBp, as in (|TBp it is ^p{p)/p^ 
that appears in the integration from zero. 

Equation (0) is similar to 't Hooft equation for a massive single flavor 

large Nc QCDg, with = It differs from 't Hooft 's equation by having 

two additional terms (two last terms in ([l^)). It suggests that the dynamics 
which governs the lowest state of the multi-flavor model is given, approxi- 
mately, by a model of massive "glueball" with an SU (Nc) gauge interaction 
and additional terms which are proportional to Nf. 

Before we present our solution of ( [l8|) it is important to note that it is only 
an approximated solution. We neglected the 3-currents state with, a-priori, 
no justification. We shall see, however, that the restriction to the truncated 
2-currents sector is an excellent approximation for the lowest massive meson. 



4 The spectrum - numerical results 



The most convenient way to solve (p!8| ) is to expand (p{k) in the following 
basis (see, however |]l5l, a different interesting choice of basis) 



^{k)=Y.Mk-7:)m-k)) 
i=0 ^ 



l3+i 



(21) 



The value of j3 chosen such that the Hamiltonian will not be singular near 
k —>■ (or k l)[|l|,|T3]. This consideration leads to the following equation 



■N, ' /3 + 1 



+ pTl cot P-K = 0. 



(22) 



This comes from eq.(|^). Had we started with (0), it would have been — /3 
replacing /? in (^5|), and constrained to /? larger than 1. Upon truncating 
the infinite sum in (pT]) to a finite sum, the eigenvalue problem reduces to a 
diagonalization of a matrix. So, the problem can be reformulated as follows 



HA 



with 



dk{k 



{k{l-k)) 



(23) 
(24) 
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and 



H., = - 1) f dk{k - \f (Ml - k)f^^^-'^^-' 



(i-M 

(k-p) 



Hence 



and 



_ B{2p + i + j + 2,2P + i+j + 2) 

2{2i3 + t+j + l) ' ^^^^ 



^ _ .iV/ i?(2/j + ^+j + l,2/? + z + j + l) 
^iV, 2(2/5 + 

B{2p + i+j + l,2p + i+j + l) 

2(2/3 + 2 + j)(/5 + J + 1) 
(/? + (/? + j)B{f3 + t,(3 + i)B{(3 + j,(3 + j) 

8(2/5 + z + j)(2/3 + ^ + j + l) 
where B{x,y) is the Beta function 



(27) 



In practice, the process converges rapidly and a 5 x 5 matrix yields the 
'continuum' results. 

The lowest eigenvalues of (0) as a function of the ratio ^ are listed in 
table 1 below (see also figure 1). Note that by /5 = 0, Nf/N^ = we mean 
the limit /3 0, Nf/N^ 0. 
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B 


Nf/N 




0.0000 





5.88 


0.0573 


0.2 


6.91 


0.1088 


0.4 


7.91 


0.1552 


0.6 


8.91 


0.1978 


0.8 


9.89 


0.2366 


1.0 


10.86 


0.2725 


1.2 


11.83 


0.3050 


1.4 


12.77 


0.3360 


1.6 


13.73 


0.3645 


1.8 


14.67 



Table 1: The mass of the lowest massive meson, in units of as a function 
of Nf/N, and /?. 



These values are in excellent agreement with recent DLCQ calculations. 

[[I2l and especially [T^ for a recent work. The typical 



For comparison see |ll 



error is less than 0.1% ! 

An interesting observation is that the eigenvalues depends linearly on Nf 
(see figure 1). The dependence is 



M^ = ^(5.88 + 5^). 

vr ^ N/ 



(29) 



We do not have a good understanding of this observation. It is not clear why 
the lowest eigenvalue sits on a straight line. It is not clear even why, as a 
eigenvalue equation (|l^) exhibits such a behavior. 

In the following sections we will consider some special cases. 



4.1 Nf = l, The 't Hooft Model 

The limit N, ^ oo with e^N, fixed and Nf <^ N, corresponds to the well 
known 't Hooft model. In this limit QCDg was solved exactly long time ago 
by 't Hooft [|1|. Let us see how our approach coincides with the fermionic 
basis in this case. In the limit Nf <^ N, we can neglect terms which are 
proportional to Nf. Equation (|T^) (or (PPD) takes the following form 
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which is just 't Hooft equation for the massless case. Note that (^) is exact, 

_i 

since in the small Nf limit the 3-currents state is suppressed by Nc ^ with 
respect to the 2-currents state and therefore we can neglect it. 

Since the wavefunction (p{k) is anti-symmetric, we will recover only the 
odd states in the spectrum of QCDg (the even states can be recovered by 
considering other sectors of the Hilbert space which decouple from the 2- 
currents state). 

Though equation (^) is formally the same as 't Hooft equation, the 
interpretation of (p{k) should be different. (p{k) is the integral of the function 
which corresponds to 2-currents state, namely to a mixture of 4-fermions 
and 2-fermions. What is the relation between the states that we find here 
and the mesons in 't Hooft 's model ? 

In order to answer this question let us expand the currents in terms of 
fermions. It is useful to denote the current in double index notation 

r{k) ^ J}{k) = dq i^¥{q)^,{k - g) - l-6^^^\q)^ ,{k - g)) (31) 

We do not bother about normal ordering, as no problem for k non zero, and 
we have to treat the k = part in a limiting way. The state | $) can be 
written as 

m = ^ / dk mJ]i-k^Jiik - 1) |0) = (32) 
dk^k) j jqdv (^¥{-q)^,{-k + q) - —5]^\-q)^u{-k + q)) X 
{^^{-p)^!,{k + P-1)- ^5l^>\-p)^k{k + p - 1)) |0). 

Note that the above expression (^2|) contains annihilation and creation fermionic 
operators. Written in terms of creation operators only, ( ^2|) reads 

|<l>) = 4r t dk t dq C ^ dp ^k)¥{-q)^>A-k + q)¥{-p)^<i{k + p-l) |0) 
2Nc Jo Jo Jo 

-7777^ f'dk t dq dp ^k)¥{-q)^,{-k + q)^^ {-p)^^{k + p - 1) |0) 
2i\^ Jo Jo Jo 

-(1 - ^) I' dk dq m¥{-q)^M - 1) |0) (33) 
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The last term in ( p3D corresponds to a meson. It can be written also as 
dq [' dk ^k)¥{-q)^i{q - 1) |0) = - C dq ^{q)¥{-q)^i{q - 1) |0) 

JO Jq Jo 

(34) 

which is exactly the 't Hooft meson. We conclude that the 2 currents state 
has an overlap with the 't Hooft meson and this is why (|18|) reproduces 
exactly the (odd part of the) spectrum of the 't Hooft model. 

4.2 Large Nf > limit 

In the limit Nf ^ N^, with e'^Nf fixed, the truncation to 2-currents state 

should again predict exact results. The reason is that the 3-currents state is 
_ i_ 

suppressed by A^^ ^ with respect to the 2-currents state. 
In this limit eq. (|T7|) takes the form 

- ^ + AV (35) 



vr \k 1 — k 

It describes a continuum of states with masses above 2m, where = 
The interpretation is clear: in this limit the spectrum of the theory 
reduces to a single non-interacting meson (or "currentball" ) with mass m. 



This phenomena was already observed in Q and in [|T7| by using a different 
approach. 

4.3 Nf = Nc, The Adjoint Fermions Model 

The case Nf = Nc is the most interesting one. It was shown that the massive 
spectrum of this model is equivalent to the massive spectrum of a model 
with a single adjoint fermion, due to 'universality' [§]. Since this model is 
not exactly solvable, it is interesting to see how our approach reproduces, 
almost accurately, previous numerical results. 

The mass of the lowest massive meson, predicted by (0), is = 10.86 x 
For comparison, the recent values reported in the literature are = 
10.8 and = 10.84 H, in units of 



This agreement is very surprising. In the regime Nf ~ N^ the 3-currents 
state is not suppressed by factors of color or flavor with respect to to the 2- 
currents state. Why, thus, is our approach so successful ? The reason seems 
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to be that as in the fermionic basis|^, the lowest massive state is an almost 
pure 2 currents state. However, the present approach is much more successful 
than the fermionic basis, where the prediction for the mass of the lowest 
massive boson of the adjoint model is twice as much as the lowest massive 
boson of 't Hooft model. It seems that the 'correct' underlying degrees of 
freedom are currents and not fermions, as predicted by the authors of P]. 

5 Summary 

In this work we used a description of massless QCD2 in terms of currents. 
With this basis we wrote down a 't Hooft like equation (|T7p for the wave 
function of the two currents states. 

The equation interpolates smoothly between the description of a single 
flavor model with large Nc ('t Hooft model), the adjoint fermions model 
Nf = Nc and the large Nf model. The equation is derived by using an a- 
priori unjustified suppression of the three currents coupling. Nevertheless, we 
observe an excellent agreement with the DLCQ results for the first excited 
state. For higher excited states the agreement deteriorates and it is of the 
order of 20%. 

The accuracy of the results for the first excited state, which implies that 
for this state the truncation of the "pair creation terms" is harmless, deserves 
further investigation. 
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Figure 1: The mass of the lowest massive meson, in units of ^-^^ as a function 
of Nf/N,. 
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